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We explore the ability of a qubit probe to characterize unknown parameters of its environment. 

By resorting to quantum estimation theory, we analytically find the ultimate bound on the precision 
of estimating key parameters of a broad class of ubiquitous environmental noises (“baths”) which the 
qubit may probe. These include the probe-bath coupling strength, the correlation time of generic 
bath spectra, the power laws governing these spectra, as well as their dephasing times T 2 . Our 
central result is that by optimizing the dynamical control on the probe under realistic constraints 
one may attain the maximal accuracy bound on the estimation of these parameters by the least 
number of measurements possible. Applications of this protocol that combines dynamical control 
and estimation theory tools to quantum sensing are illustrated for a nitrogen-vacancy center in 
diamond used as a probe. 


I. INTRODUCTION 

Controlled spin-i particles (qubits) are sensitive 
probes of the structure and properties of highly complex 
molecular, atomic or solid-state quantum systems. Novel 
quantum technologies requiring such high sensitivity at 
the nanoscale are based on qubit probes serving as sen¬ 
sors [1-8] or monitors of biological or chemical processes 
[9-11]. Here we focus on the extraction of information 
characterizing the environment of a qubit probe, by mon¬ 
itoring the decoherence process that the qubit undergoes 
[12]. Dynamical control, originally conceived as a tool for 
reducing decoherence effects [13-19], has also been shown 
to be a valuable source of information on environmental 
noises [10, 20-25]. This information is revealed by the de¬ 
pendence of the decoherence rate of the qubit probe on 
a control-field parameter, owing to the fact that, under 
weak coupling of the qubit and the environment (“bath”), 
this rate is universally expressed by the overlap of the 
environmental noise spectrum and a spectral filter func¬ 
tion that is solely determined by the dynamical control 
[15-17, 26]. The filter function can then be designed by 
varying the control field to scan the noise spectrum. The 
information obtained from this procedure, dubbed “noise 
spectroscopy” [22, 23], is not only essential for designing 
the most effective (optimal) dynamical protection from 
decoherence caused by a given environment [26-33] of 
quantum-information processing [34-37], quantum-state 
transfer [38, 39], Hamiltonian engineering [24, 38, 40] and 
quantum-state storage [41, 42]. It may also become a 
means of understanding physical or chemical processes 
[9-11, 24] by analyzing their noise fluctuations in mag¬ 
netic resonance spectroscopy and imaging with nanoscale 
resolution [6-8, 11, 43-45]. 

In order to further advance this promising noise spec¬ 
troscopy and broaden its applicability, it is imperative 
to find the best general strategy for extracting the envi¬ 
ronmental (“bath-induced”) noise-spectrum information 
from the qubit-probe decoherence. The strategy we 
adopt aims at minimizing the error in estimating key pa¬ 
rameters of the noise (bath) spectrum by measurements 
performed on the qubit probe (Sec. H). The minimal er¬ 


ror is determined by the maximal quantum Fisher infor¬ 
mation (QFI) [46-49] gathered by measurements in the 
optimal basis (Sec. HI). We here hnd (Sec. HI A) the ul¬ 
timate error bounds for unbiased estimators when the 
qubit probe, under arbitrary control, undergoes pure de¬ 
phasing in the probe-bath weak-coupling regime. Moti¬ 
vated by practical experimental considerations and con¬ 
straints, these ultimate bounds attain the best estima¬ 
tion precision by the least number of measurements pos¬ 
sible. We achieve these goals by replacing the free- 
evolution (-induction) decay (FID) of the qubit state 
prior to each measurement by dynamically-controlled evo¬ 
lution designed to ensure the convergence to the ultimate 
error bound for the particular bath-spectrum parameter 
to be estimated (Sec. HIB). For each such parameter, an 
appropriate filter-function must be generated by dynam¬ 
ical control (Sec. IV A) [17, 26]. We here focus on deter¬ 
mining the general conditions that have to be satisfied 
for designing the filter function to attain the ultimate 
bounds. 

The first step in the proposed strategy is the estimation 
of the coupling strength g of the bath to the probe, that 
can be interpreted as the noise variance. The proposed 
appropriate filter function is generated by projections 
onto an eigenstate of the Ux probe-operator at a rate 
that conforms to the quantum Zeno regime [15, 24, 50] 
(Sec. IVB). This procedure, which does not require prior 
knowledge of the bath spectral lineshape, has been experi¬ 
mentally exploited to determine the coupling strengths of 
complex spin-networks without maximizing or optimizing 
the information obtained [24]. 

Once the coupling strength g is optimally estimated, 
the bath spectra have to be characterized by their nor¬ 
malized lineshapes. These spectra crucially depend on 
the bath correlation time Tq often unknown to us: it is 
typically the inverse of ujc, which is the width or the cutoff 
of the bath spectrum. We show that convergence to the 
lowest error bound on Tc is achievable through filter func¬ 
tions generated by common types of coherent dynamical- 
control sequences (Sec. IV C). By contrast, free-induction 
decay (FID) of the qubit coherence {(Jxit)) may preclude 
such convergence. 
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We further show for a family of generic baths that op¬ 
timized convergence requires a filter function that only 
samples (overlaps) a power-law region of the bath spec¬ 
trum (Sec. Ill B). Such spectral features characterize om¬ 
nipresent baths: sub-Ohmic, Ohmic and super-Ohmic 
baths whose spectra obey a power-law at low frequen¬ 
cies, as well as noise spectra of generalized Ornstein- 
Uhlenbeck processes characterized by a power-law tail at 
high frequencies. These types of bath spectra are ubiq¬ 
uitous in solid-state, liquid or gas phases [11, 23, 51- 
53] where they represent collisional or diffusion processes 
[54, 55]. Other environmental parameters, such as spec¬ 
tral power-law exponents, the T 2 decoherence time and 
diffusion coefficients, are shown to obey bounds analo¬ 
gous to those of g and Tc (Sec. IVD). 

Finally, we demonstrate the practical feasibility of ex¬ 
periments that may attain the ultimate precision bounds 
by resorting to a real-time adaptive estimation proto¬ 
col based on a Bayesian estimator and an online exper¬ 
imental learning design [56-59] (Sec.V). We resort to 
this protocol to illustrate the ability to achieve the pre¬ 
dicted analytical bounds in an efficient way under exper¬ 
imentally relevant conditions, e.g. for nitrogen-vacancy 
center (NVC) probes in diamond [1-8, 45, 51]. Thus, 
the present analytical theory, supported by adaptive- 
estimation simulations, suggests that the proposed syn¬ 
thesis of optimally-controlled noise spectroscopy and es¬ 
timation theory can become a powerful, broadly applica¬ 
ble, diagnostic tool (Sec. VI). 


II. DEPENDENCE OF THE QUBIT-PROBE 
DEPHASING ON AN UNKNOWN 
ENVIRONMENTAL PARAMETER 

We consider a qubit-probe that experiences proper¬ 
dephasing by the environment (bath) under the action 
of the system-bath interaction Hamiltonian 

Hsb = gcTzB, ( 1 ) 

where g is the probe-bath interaction strength, az is the 
appropriate Pauli operator for the probe and B is the 
bath operator (Fig. la). To obtain maximal information 
on the environment, a convenient initial probe-state is the 
symmetric superposition of the qubit-up/-down states in 
the az basis, 

^(it)+ !;)) = !+), (2) 

and the optimal observable for the probe-state measure¬ 
ment is ax [49] (Appendix B). Our goal is to determine 
the dynamical control efficacy for estimating, one by one, 
the unknown parameters such as g and Tc that charac¬ 
terize the coupling bath spectrum. 

We shall denote the particular bath parameter by 
xb, and, accordingly, the bath coupling spectrum by 
G{xb,ui). Under proper dephasing, 

{ax{xB,t)) = Tr{ps{xB,t)ax) = (3) 


where fr{xB,t) is the attenuation factor due to dephas¬ 
ing. In the probe-bath weak-coupling regime known as 
the Born approximation for the bath, wherein the qubit- 
probe negligibly influences the environment [60] , it obeys 
the universal formula (Appendix A) [15-17, 26] 

/ OO 

duiFt{uj)G{xB,wi), (4) 

-00 

where Ft(u}) is a filter function which explicitly depends 
upon the dynamical control of the qubit-probe during 
time t. 

This universal formula is exact for Gaussian noise. Yet 
it applies to any noise under the weak-coupling assump¬ 
tion in Eq. (4). In our bath-optimized control theory 
[15-17, 26, 27], the filter function is obtained to ensure 
optimal control for any given bath and task at hand. This 
theory does not impose any requirement on the temporal 
shape of dynamical control that defines the hlter func¬ 
tion: the control may be continuous or pulsed, coherent 
or projective, in contrast to the dynamical-decoupling 
formulation [14, 18, 29]. We stress that the optimal filter 
function obtained by this bath-optimized control theory 
does not require high-order corrections as opposed to dy¬ 
namical decoupling control [26, 27]. 

The information about the unknown bath parameter 
Xb is encoded in the protocol dehned by Eqs. (l)-(4) 
(Eig. la) by the probabilities p of Ending the qubit in 
the |-|-) (symmetric) or |—) (antisymmetric) state when 
measuring ax- These probabilities obey 

p(±|xB,t) = i(l±e-‘^("-’‘)). (5) 


HI. OPTIMAL ESTIMATION UNDER 
DYNAMICAL CONTROL 


The minimum achievable relative error of the unbiased 
estimation of a single unknown parameter xs is deter¬ 
mined by the quantum Cramer-Rao bound to be 


e{xB,t) = ^^> - ^ ^=, ( 6 ) 

Xb XB\/NmFQ{xB,t) 


where Nm is the number of measurements and Fq{xb, t) 
is the quantum Fisher information (QFI) that quantifies 
the maximum amount of information on xb that can be 
extracted from a given state [46]. Therefore, we set out 
to maximize the information 


Fq{xb 1 ^opt ) = maxJ'e(xB,t), (7) 

by choosing the optimal time topt to perform the mea¬ 
surement and an appropriate dynamical-control scheme 
(prior to the measurement), to obtain the minimal attain¬ 
able (relative) error in the estimation of xb for a given 
bath spectrum. For a qubit-probe obeying Eq. (5), the 
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Figure 1: Scheme of noise-spectra parameter estimation, (a) Scheme of the probing process and its dynamical control: Dy¬ 
namically controlled qubit-probe undergoes dephasing by the environment (bath), (b) An optimal qubit-probe filter function 
generated by frequent projections is spectrally flat to the extent that it conforms to be quantum Zeno regime and allows 
the optimal determination of the coupling strength g regardless of the shape of S{uj). (c) A filter function of an optimally 
controlled qubit-probe, for obtaining the ultimate bound on the correlation time (tc) estimation. Such a filter must 

overlap only with the power-law tail of the Ornstein-Uhlenbeck noise spectrum Sg{Tc,to), or with oj > 0 but not with tUc for 
the (super) Ohmic spectrum Sa{rc,ijj). The free-evolution (FID) filter function Fjf'~‘^'^{ijj) does not fulfill the requirements for 
achieving the bound, since it is centered at cu = 0. 


QFI is gathered by measurements in the optimal basis 
(Appendix B) 


an exponent a, where the derivative of J7 ’(xb, t) is tightly 
bounded by its value in the power-law region, as 




1 — ^ dxB ) 


( 8 ) 


dJ{XB,t) 

dXB 


< 


aJ{xB,t) 

xb 


( 10 ) 


From this expression, the QFI maximum is obtained by 
finding the optimal tradeoff between the signal-amplitude 
contrast (the magnitude of ) and the sensitivity of 

the signal-attenuation to the parameter Xb (the deriva- 
tive ■§^)- Obviously, neither should be too small if Eq. 
(8) is to be maximized. Their optimal tradeoff determines 
topt for a given control scheme (Fig. 2a). To evaluate the 
efficacy of different dynamical controls in attaining the 
highest accuracy in xb, we then use the error (6) at topt, 
i.e. s{xB,topt), as the figure of merit. 

We focus here on the ultimate error bounds for the es¬ 
timation precision per measurement for the qubit-probe 
that undergoes pure dephasing within the weak-coupling 
regime. These bounds are of practical importance in typ¬ 
ical experimental situations where the initialization and 
readout times, normally determined by Ti ^ T 2 , con¬ 
strain the time interval between measurements. 

Alternatively, in instances where Ti is not the domi¬ 
nant constraint, one may be interested in the ultimate 
precision bound attainable during a given interrogation 
time that extends over many measurements, T = Nmt. 
In such cases the Fisher information per unit time has to 
be maximized: 

J^Q{xB,topt)ltopt = max{TQ{xB,t)/t). (9) 


A. Ultimate estimation bounds 

As discussed below, a broad class of practically relevant 
environmental noise processes cause the attenuation fac¬ 
tor J{xB,t) to have a power-law dependence in xb with 


The equality strictly holds when J{xB,t) is an homoge¬ 
neous function of degree a for the parameter x b ■ 

The bound (10) leads to a tight lower bound on the rela¬ 
tive error in (6), which holds for either free decay (FID) 
or arbitrary dynamical control of the probe coherence 
(Appendix C) 


e{xB,t) > 


Vl — 

ay/Nm Jo e-'Fo 


£0 


( 11 ) 


Here the value of J = Jo that minimizes this bound, 
Jo = 1-h ^W(—2e~^) ~ 0.8, W{z) being the Lambert 
function, yields the equality in (11), with eo ~ 2.48. We 
have thus obtained a universal, ultimate error hound for 
the unbiased estimation of xb in a broad class of noise 
spectra probed by a qubit-probe undergoing dephasing 
within the probe-bath weak-coupling regime. 


B. Key parameter estimation 

This general bound applies to the estimation of any pa¬ 
rameters that satisfy J {xb , t) oc . We here discuss in 
detail some examples of estimating key parameters that 
belong to this class in typical scenarios of practical im¬ 
portance: 

(i) The effective coupling strength g of the qubit-probe 
with the dephasing bath is required for defining 

G{xB=g,uj) = g'^Siuj), ( 12 ) 

S{lo) being the normalized spectral density of the bath 
autocorrelation function, so that G{g,uj) and J{g,t) are 
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Figure 2: Optimization for achieving the ultimate precision bound, (a) Magnetization {a^) as a function of the measurement 
time t of a qubit-probe experiencing dephasing due to Ornstein-Uhlenbeck process, with Lorentzian spectrum G^= 2 (tc = 10, w) 
with g = 1 under a CPMG control sequence with N = 8 pulses using the exact analytical expression derived from Eq. (4) [11]. 
The optimal time that determines the best tradeoff between the signal contrast (almost halfway between 0 and 1) and the 
highest sensitivity of the decay rate to Tc, provides the most accurate estimation of Tc- (b) Minimal relative error 
per measurement {Nm = 1) of the estimation of Tc as a function of gr^ for a Lorentzian spectrum G^= 2 (tc = 10, tj) (dashed 
lines) and a super-ohmic spectrum Gb^ 2 (tc = 10, a;) (solid lines) calculated from the integral of Eq. (4), both with g = 1. The 
minimal relative error under optimal control (CPMG/CW with N pulses/cycles satisfying the conditions given in the main 
text, green lines) achieves the ultimate bounds (dashed-) and (solid-) green lines for s = 2 and /3 = 2 respectively). 

It can improve the best precision obtainable under free-evolution (blue lines) by several orders of magnitude. 


homogeneous functions of degree a = 2 \n g. 

(iia) The correlation time Tc of noise fluctuations is 
required for describing generalized Ornstein-Uhlenbeck 
processes, with normalized spectral densities 

Sf}{xB=Tc,uj) = Ap -(13) 

1 + UJ^Tc 

where /3 > 2 is an even integer and Ap = is the 

normalization factor. 

(iib) The correlation time Tc = that is the inverse 
of the cutoff frequency Wc, in generalized Ohmic spectra 

Ss{xb=Tc, w) = (s -I- l)w“<^'’+^)u;®0(a;)6(a; - Wc), (14) 


where ©(w) is the step function, s = 1 stands for a Ohmic 
spectrum, 0 < s < 1 for a sub-Ohmic spectrum, and 
s > 1 for its super-Ohmic counterpart. 


Both spectral densities (13) and (14) satisfy 


dS 

OTc 


< 


and consequently the bound (10). Furthermore, 
both spectra attain the ultimate bound (11) at fre¬ 
quency ranges where they have power-law dependence 
(Appendix C): spectral density (13) at high frequen¬ 
cies, where it becomes a homogeneous function of degree 
a = /3 — 1, and spectral density (14) at low frequen¬ 
cies, when we restrict ourselves to to < lOc to avoid the 
cutoff effects, thus rendering the spectral density a ho¬ 
mogeneous function of degree a = s -I- 1. 

(iii) The dephasing time T 2 in the attenuation exponent 
which is a homogeneous function of degree a, J(xb = 
= (^"^ 2 ) , may be estimated down to the ultimate 


bound. 


IV. DYNAMICAL CONTROL STRATEGIES 
FOR ACHIEVING THE ULTIMATE BOUND 

A. Achieving the ultimate-bound by dynamical 
control 

It follows from the discussion above that the ultimate 
error bound, Eq. (11), in the estimation oi xb = g or Tc, 
can be attained provided that the dynamical control on 
the probe generates a filter Ft{uj) that extends only over 
the frequency band where the noise spectrum behaves 
as a power-law in and the equality in Eq. (10) is 
fulfilled. Then, upon adjusting the total control time 
such that J[xbA°^*) = the equalities in Eq. (11) 
are also fulfilled. 

We note that these are general conditions for attaining 
the ultimate error bound per measurement, but they do 
not invoke the optimization of the filter function under 
specific constraints that may be imposed in a given ex¬ 
perimental setup. In such setups further optimization is 
required to approach as best we can the ultimate error 
bounds. 

B. Estimating the probe-bath interaction strength g 

If 5 ( 0 ;) is unknown apart from a crude estimate of its 
overall width, dynamical control is needed for estimating 
g down to the ultimate precision bound. The appropriate 
control is such that the filter function Ft{uj) = Ft is flat 
in the domain of S'(a;), leading to an attenuation factor 
ff = Ftg"^ that is independent of the noise spectrum. For 
free evolution this limit holds when the interval is much 
shorter than the correlation time Tc of the environment, 
so that the qubit-probe evolves freely (Appendix D). This 
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condition can only be fulfilled for large enough Tc- 

To overcome this limitation, we may instead realize 
such a filter at times larger than Tc via repeated strobo¬ 
scopic projections of the qubit-probe in the basis of |±), 
i.e. by quantum non-demolition (QND) measurements 
at a rate that is high enough to conform to the quan¬ 
tum Zeno regime [50]. The filter function then describes 
spectral broadening of the |±) eigenvalues far beyond the 
width of S{uj) (Fig. lb). The attenuation factor is then 
[15, 17, 24] (Appendix D) 


Jzenoig^t) = if (15) 

obtained for N QND measurements during the total con¬ 
trol time t. The advantage of this Zeno regime is that 
it requires ^ r^, rather than the total time t, to be 

less than Tc- Since the outcomes of these measurements 
need not be read out, but rather used to guide the evo¬ 
lution, they may be emulated by impulsive noise-induced 
dephasing of the qubit probe that has the same effect as a 
projection on the probe evolution [17, 24, 61]. This QZE 
regime has been already exploited experimentally to de¬ 
termine the coupling strengths of complex spin-networks 
[24]. 

Under the quantum-Zeno condition of (15), the error 
estimation bound (11) is achieved when the total control 
time (after which ax is measured and read out - see (5)) 
is chosen to have the optimal value 

, provided that N > (16) 

This equation (further discussed in Appendix D) consti¬ 
tutes the main condition on an optimal filter designed to 
estimate g. 

C. Estimating the correlation time Tc 

To achieve the ultimate precision bound (11) for the es¬ 
timation of Tc of noise spectra (13)-(14), the control on 
the probe should generate a hlter Fi(w) that only overlaps 
with the power-law portion of 5^(8) (tc, w) oc ^ 

as shown in Fig. Ic. By contrast, FID of the probe co¬ 
herence generates a filter (sine) function centered at zero 
frequency, thus preventing the hound in Eq. (11) from 
being reached. 

While various controls may allow the best estimation ac¬ 
cording to Eq. (11), we here analytically study the con¬ 
ditions for achieving the ultimate bound under standard 
CPMG sequences of tt pulses [62], or under continuous- 
wave driving (CW). In control sequences of equidistant 
pulses N :s> I, the filter function Ft{oj) converges to a 
sum of delta functions (narrowband filters) centered at 
the harmonics of the inverse CPMG time period, while 
for CW there is a single frequency component [17, 23, 30] 
(Appendix D). In the following we use these filter func¬ 
tions, to analytically infer the required total control time 


t and N that allow the bound to be attained for the two 
classes of power-law spectra in Eqs. (13) and (14) (see 
Appendix D for details): 

(i) For generalized Ornstein-Uhlenbeck spectra Sp{uj), 
only high frequencies must be probed by the dynam¬ 
ical control filter. To this end, the intervals between 
pulses or refocusing periods must obey ^ Tc- Then, 

flp (X Tc satisfies the equality in (11) and the min¬ 

imal relative error attains the ultimate bound, provided 
that the total control time is chosen to have the optimal 
value 


.opt 

*/3 


= Tc 



9 9 ' 

cpg Tf 


(17) 


where cp is a constant depending on the control sequence. 
The bound can be attained only if N is sufficiently large 

jopt 

to satisfy t°^ , Tc ^ and only overlaps with the 
power-law tail, as shown in Fig. Ic. By contrast, FID 
fails this condition, since its filter mainly overlaps with 
a; ~ 0 (Fig. Ic), causing a larger error in the estimation, 
as shown in Fig. 2b. 

(a) For generalized Ohmic spectra Ss, the filter should 
only overlap with the noise spectrum at a frequency lower 
than the cutoff, Q < ut < Uc, avoiding any overlap at Wc. 
Then, Js oc satisfies the equality in (11) and the 
minimal relative error attains the ultimate-bound when 
the measurement of ax is performed at the optimal total 
control time 


.Jcsg'^T^N'^ 

V Jo 


(18) 


topt 

Here > Tc, and N ^ 1 ensure the filter is narrowband 
with negligible tail at oje, as opposed to its free-evolution 
(FID) counterpart that causes a much larger estimation 
error (see Fig. 2b). 

Equations (17) and (18) thus constitute our main con¬ 
dition on the optimal filter design for estimating Tc- 


D. Power-law estimation 

One can optimally estimate the exponent /3 or s that 
governs the bath spectrum by maximizing the QFI (8) for 
the estimation ofa;B =7 = /3-|-l = l — s (Appendix D). 
This maximization leads to the ultimate precision bound 

(- w 1 \ Vl — e~^Fi 

^ “ 7V^-^C(7U) “ VK^Ji\ln{Ji)\ e-^^^7^ 

(19) 

with Ji = e Ri 0.106 and ei ~ 2.04. The bound 
is achieved when the qubit only probes the power-law 
regime of noise spectra (13)-(14). Provided that the 
GPMG or CW dynamical control probes the power law 
region, the bound is attained when the qubit-probe is 
measured at 


topt — T2 \fjl, 


( 20 ) 
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where and T 2 is the dephasing time 

(that depends on the applied control). 

V. REAL-TIME ESTIMATION PROTOCOL 

The predicted optimal time for performing a mea¬ 
surement following the dynamical control or FID of the 
quantum-probe coherence, explicitly depends on the un¬ 
known parameter xb to be estimated (see Eqs. (16), (17), 
(18) and (20)). In practice, one may bypass this difficulty 
by estimating xb and simultaneously finding the optimal 
time to monitor the probe via an efficient real-time es¬ 
timation protocol [56, 57, 59]. We illustrate here the 
implementation of this protocol and show that it attains 
the predicted ultimate bound (11) if the qubit-probe is 
optimally controlled. 

The protocol is based on a Bayesian estimator and an 
online experimental learning design that maximizes the 
information gain, similar to the one that was recently im¬ 
plemented for improving the suppression of decoherence 
in quantum dots [58]. Its stages are as follows: 

(i) We initially guess a probability distribution p(xb) for 
the unknown parameter xb that represents our apriori 
knowledge of the parameter to be estimated. Physically, 
we should have > 0, and therefore one can assume a 
flat distribution oi p{xb), xb €1 {x‘]f™, x‘]f°'^] that con¬ 
tains the true value 

(ii) We then determine the best time tm to perform a 
measurement for maximizing the information gain about 
Xb which is defined by the information entropy Ultm) = 
max{U{t)'\ (Appendix E), which depends on p(xb) and 

the likelihood function p{d\xB,t) of Eq. (5) that deter¬ 
mines the conditional probability to obtain the possible 
outcome data of the measurement d = {-I-, —}. 

(Hi) We next perform a measurement of ax of the qubit- 
probe state, initialized in the state |-|-), at the optimal 
time tm, obtaining the outcome data d with probability 
p(d|x^“®,t). According to the obtained data, we update 
our knowledge of by the Bayesian Rule 

p[d\tm) 

where p{d\t) is a normalization factor for integration over 
Xb- 

The estimation oixB improves upon iteratively repeating 
this three-stage process Nm times, Nm standing for the 
number of measurements. When an adequate control is 
chosen, the probability distribution p{x b ) converges to a 
narrow peak around By contrast, the convergence 

under free-evolution can be very poor. 

Figure 3 presents a simulated experiment of this itera¬ 
tive process for the estimation of g and Tc of the envi¬ 
ronmental noise by a NVC spin-probe in two types of 
diamond samples whose intrinsic environmental charac¬ 
teristics require dynamical control to attain the ultimate 
bound (11): 


(a) We simulate the estimation of Tc for the NVC-probe 
evolving freely (undergoing FID), as compared to it being 
nearly optimally controlled by a CPMG sequence with 
N = 8 subject to an Ornstein-Uhlenbeck process charac¬ 
terized by G/3=2('rc'^“'^ = 10/rs, w) with g = 1 MHz (con¬ 
sistently with the HPHT diamond sample data of Ref. 
[51]). Our simulated best-measurement timings tm con¬ 
verge to the theoretically predicted t°P* that maximizes 
the QFI (8), as exemplified in Fig. 3a for CPMG control, 
where t°P* k- 18.25^s is on time scales compatible with 
the accessible experimental times [51]. Concurrently, the 
minimum relative error converges to e{Tc,topt) predicted 
from the Cramer-Rao bound (Fig. 3b). 

(b) In Fig. 3c we simulate the estimation of g under 
similar conditions by a sequence of iV = 500 QND mea¬ 
surements due to induced dephasings whose optimal in¬ 
tervals are « 1.9^s and compare these results to the 
FID results. In this case G = 0.03MHz,a;) with 
Tc = lO^s is consistent with the spectral density deter¬ 
mined in Ref. [51] for the ^^C diamond sample. This 
estimation procedure is performed without assuming the 
lineshape of the spectral density. It is seen from Fig. 3c 
that the precision error eventually saturates for the free 
evolution (FID) at a precision that is limited by our lack 
of knowledge of the lineshape. By contrast, this satu¬ 
ration can be overcome in the QZE regime achieved by 
projective measurements. 

The foregoing simulations have thus confirmed, under ex¬ 
perimentally realistic conditions, our second major ana¬ 
lytical result, whereby the ultimate theoretical bound (11) 
is indeed attainable under optimal control, but hardly ever 
under free-evolution. 


VI. DISCUSSION 

We have demonstrated that dynamical control of a quan¬ 
tum probe not only dramatically improves the quan¬ 
tum estimation of environmental parameters compared 
to the free-evolution (induction) decay (FID) of its co¬ 
herence: it may be imperative to use such control, since 
FID may preclude their correct estimation. In particu¬ 
lar, for generic noise spectra, as in generalized Ornstein- 
Uhlenbeck and Ohmic processes, the ultimate analytical 
bounds for the coupling strength and correlation-time es¬ 
timation precision derived here can be achieved by opti¬ 
mizing the dynamical control on the probe. The opti¬ 
mal controls suitable for the estimation of g, Tq (Fig. 
1, 3) or the power-law exponent are generally different, 
but the protocol is similar. Specifically, we demonstrated 
that the ultimate estimation bound for g can be attained 
in the quantum Zeno regime without prior knowledge 
of the environmental spectral-density shape. Once this 
probe-environment coupling is known, the environmental 
correlation time can be reliably estimated by standard 
dynamical-control sequences. 

Experimental conditions typically require a minimized 
average error within the coherence time interval shown 
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Figure 3: Simulation of a experimental real-time adaptive estimation protocol for realistic conditions with a NVC spin-probe. 
(a,b) Convergence of the real-time adaptive estimation protocol to the theoretically predicted values for estimating Free- 
evolution of the probe (blue circles) is contrasted with that of a dynamically controlled probe under CPMG (green square) 
sequence with = 8 in the presence of an Ornstein-Uhlenbeck process with Lorentzian spectrum G/3=2(Tc’'“'^=10^ts, ca), with 
g = 1 MHz consistently with the spectral density of a HPHT diamond sample determined in Ref. [.51]. The simulated curves 
derived from exact analytical results of Eq. (4) [11], were averaged over 600 realizations. In (a) the optimal measurement time 
tm as a function of Nm converges to the value for the CPMG case. Similar curves converging to the corresponding 
are observed for other controls and free evolution. In (b) the minimal relative error converge to the (Cramer-Rao) 

bound. Under free evolution the regime where e oc is attained for Nm ^ 100. The ultimate bound (dashed line, 

a = j3 — 1) is only attained by optimal control, (c) Convergence to the minimal relative error e{g,topt) to the (Cramer-Rao) 
bound for estimating g hy N = 500 consecutive projective measurements in the Zeno regime (green triangle) compared to the 
estimation under free evolution (blue circle). In this case G/3=2((?=0.03MHz, w), with Tc = lO/rs, consistently with the spectral 
density of a ^^C diamond sample determined in Ref. [51]. Here too the ultimate bound dashed line, a = 2) is only 

attained by optimal control, (d) Proposed scheme for using a NVC as a qubit-probe for its environment. The ms — 0 ([0)) 
state is fully populated by laser irradiation (dashed curly arrow). Microwave (MW) pulses are selectively applied between the 
states with = 0 and —1 ([0) and [ — 1)) to initialize the spin-probe in a [-|-) = -^([O) + [ — 1)) state, and to apply the tt pulse 
CPMG sequence for estimating t^. For estimating g, projective measurements are performed by combining MW 7 r/2 pulses 
on the 0 o — 1 transition and laser-induced relaxation between the ground and exited electronic states that conserve the spin 
components (solid curly arrows). The readout is done at the end of the sequence by detecting the laser-induced fluorescence 
signal. 


in Fig. 2b, if this interval is much shorter than the com¬ 
bined post-measurement readout and initialization time 
of the qubit. Under such conditions, an important impli¬ 
cation of the present analysis is that it yields the opti¬ 
mal total time of the pre-measurement dynamical control 
and thereby the number of pulses and the delay between 
them that need be applied to attain the best precision per 
measurement. However, in the (rarely encountered) op¬ 
posite limit of fast readout and initialization, the overall 
error rate of the measurements is to be minimized. Then, 
a similar optimization can be performed by maximizing 
the Fisher information per unit time. If additional spe¬ 
cific constraints are imposed by a particular experimental 
setup, e.g. limited power or estimation time, then further 
optimization is needed to find the optimal filter function. 
In such cases, the goal would be to approach as best we 
can the ultimate error bound under the given constraints. 


A real-time adaptive estimation protocol has illus¬ 
trated the ability to find the optimal time of dynamical 
control for achieving the predicted precision bounds, e.g., 
for a nitrogen vacancy center (NVC) in a diamond that 
acts as a qubit-probe of the bulk or surface in optically 
detected NMR and MRI. Its aim here is to determine the 
noise spectrum generated by nuclear or electronic spins 
[6-8, 44, 45], as well as other sources of dephasing. In Fig. 
3 we illustrated this protocol for estimating both g and 
Tc near an NVC used as a probe. The inferred correlation 
times of the noise fluctuations can be particularly helpful 
for studying molecular diffusion at the nanoscale, where 
the power law tails of noise spectra allow the extraction 
of the diffusion rates and restriction lengths associated 
with pore structures and thereby characterizing biologi¬ 
cal processes [11, 43]. 

In other scenarios, power-law tails can also characterize 



































charge diffusion in conducting crystals [55] or spin diffu¬ 
sion in complex spin-networks [23, 24, 54, 63]. General¬ 
ized Ohmic spectra may help understanding the function¬ 
ing of nanoscale electromechanical devices [64], as well 
as superconducting devices attached to conducting leads 

[65]. 

Several extensions of the outlined strategy will be fur¬ 
ther explored: 

(a) If uncontrollable noise sources are present, such as 
an intrinsic T 2 decoherence due to a white-noise (Marko¬ 
vian) process or due to pulse imperfections, we may have 
to resort to more elaborate controls, nicknamed by us 
“selective dynamical recoupling” (SDR), which we have 
shown to allow selective probing of the targeted noise 
source [10, 11]. SDR was already implemented to se¬ 
lectively monitor diffusion processes characterized by an 
Ornstein-Uhlenbeck spectral density, so as to determine 
the probe-environment coupling strengths [10, 66] and 
the environmental correlation time, as well as diffusion 
restriction lengths [11, 43] in systems of biological and 
chemical interest. We envisage that by incorporating 
SDR within the present optimized estimation strategy 
one may eliminate Markovian (intrinsic-T 2 and pulse- 
error) effects and allow a clean parameter estimation of 
the targeted noise spectrum. 

(b) A qubit-probe may be replaced by an entangled n- 
particle probe that may yield a lower error bound [67-69] . 
The estimation strategy will remain similar, but the dy¬ 
namical control will have to be adapted to a multipartite 
scenario, based on the approach of Refs. [26, 70, 71]. 

(c) The present strategy, whereby environmental pa¬ 
rameters are estimated one by one, provides the ultimate 
precision bounds that can be attained for estimating a 
single parameter, assuming that the remaining parame¬ 
ters are known with much better accuracy. Therefore, 
these bounds are also valid for more elaborate multipa¬ 
rameter estimation strategies based on a quantum Fisher 
information matrix [72] , although these bounds may not 
be tight in the latter case. Still, the single-parameter 
estimation expounded here may be superseded by multi¬ 
parameter optimized estimation per measurement. 

To conclude, we have analytically set the error bounds 
on environmental-parameter estimation and demon¬ 
strated the ability to extract such environmental infor¬ 
mation with maximal accuracy by the least number of 
measurements possible upon invoking dynamical control. 
This general result for a qubit-probe weakly coupled to 
its environment opens the door to the development of an 
important diagnostic tool of environmental processes by 
quantum probes. 


Appendix A: Anamnesis: Filter function derivation 
for qubit probe under dephasing (Following Refs 
[14-16,24-25]) 

The Hamiltonian of a single spin-probe interacting with 
a bath that produces pure dephasing is 

H = Hs{t)+HB+HsB, (Al) 

with 

Hs = YCT. + /(t)e-*‘"“V„ Hsb = S®B = ga,B, (A2) 

where f(t) is the dynamical control applied to the qubit, 
B and S are the operators of the bath and the system (cr^ 
for pure dephasing) respectively, and g is their coupling 
strength. 

In the interaction picture the Hamiltonian can be writ¬ 
ten as 

H^SBit) = S{t)(g,BHt), (A3) 

where 

S{t) = uUt)SUs{t), Us{t) = 

Bit) = Ul{t)BUB{t), UBit) = 

Therefore, Eq. (A3) becomes 

Hsb = 9^it)(^zBit), (A5) 

where D(f) is the dynamical control rate of the system. 

From this form one can derive the non-Markovian mas¬ 
ter equation for the density matrix of the system, psit) 
in the interaction picture, which in the Born approxima¬ 
tion one assumes a weak coupling g such that the the 
influence of the qubit-probe on the environment is small 
(usually called the weak-coupling approximation). In this 
approximation, the density matrix of the environment 
Pb is only negligibly affected by the interaction with the 
qubit-probe, and the state of the total system at time t 
is allowed to be expressed as p{t) ps{t) ^ pb [60]. The 
resulting non-Markovian master equation is then given 
by [15-17, 73, 74] 

t 

PsiTc.t) = J dt'{g^<^ixB,t - t')[Sit'), Sit)psit)] -\-h.c}, 
0 

(A6) 

where 4>(xB,f'—f”) = Tr^ {B{t' — f”)i?(0)pB(0)} are the 
bath correlation functions and xb is a parameter that 
characterize the environment. Then, the attenuation 
factor of the spin coherence (magnetization) {oxit)) = 
nieasured at time t is 

t t' 

JixB,t) = j dt' j df^g^<^ixB,t'-r)nit')n*if^), 

0 0 

(A7) 
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which can be cast in the spectral form 


/ OO 

duiFt{uj)G{xB,(^), 

-OO 


(A 8 ) 


1 f°° 

G{xB,i^) = Tr dtg^^{xB,t)F 
7_oo 

being the bath-coupling spectrum and 

/ dt'n{t')F‘^* 

Jo 


Ft{uj) = 


1 

2tt 


(A9) 


(AlO) 


is the filter function which depends on the dynamical 
control of the qubit probe. 

We stress that Q{t) can have an arbitrary temporal 
shape, as opposed to the restrictions on its {inherently 
pulsed) shape in the dynamical decoupling method (Refs. 
[12,13,17,27] in the main text). 


Appendix B: Quantum Fisher information 
concerning a single parameter of the environment 

The quantum Fisher information (QFl) concerning a sin¬ 
gle parameter xs of the environment (bath) for the qubit 
state is [46, 47, 49] 




p+ \dxB J P- xdxB/ 


, r,iP+-P-? 
^ P++P- 



d\p+) 

dXB 


(Bl) 


(P+I 


d\p-) 


dxi 


where 

p±{xBt) =p{±\xB,t) = i (^1 ± , (B2) 

b±)=^(e—“‘mill)) 

are the eigenvalues and eigenvectors of the spin-probe 
density matrix [17, 74]. A measurement is said to be 
optimal when the QFl Fq coincides with its classical 
counterpart [46, 47, 49]. This is the case here under 
pure dephasing, when the last term in (Bl) is null due to 
= 0. Then the optimal measurement is effected by 
projections onto the eigenstates |p±) of <Jx in the rotating 
frame 


\p±){p±\=\e-^^-^^^ |±) 

|±) =y=(|t)±|;)). 


Correspondingly, Eq. (Bl) becomes 


FQ{xB,t) 


e 2J(2:i3,t) /dJ{xB,t)\^ 

1 — \ dxB ) ’ 


(B3) 


(B4) 


if the initial probe-state is |-|-). An arbitrary initial 
state, (cos(0) It) + * sin(0) ID), 0 < 6 < ^, leads to 
FQ{xB,t) oc sin^(20) [49]. Therefore, the optimal initial 
state leading to the maximal QFl is obtained for 0 = j, 
|-|-) = ^ (It) + It)) , thus proving Eq. ( 6 ) in the main 
text. 


Appendix C: Derivation of the ultimate precision 
bound 


A broad class of practically relevant environmental 
noise processes cause the attenuation factor to have a 
power-law functional dependence in xs with exponent 
a, where the derivative of J{xbF) with respect to xb 
satisfies 


dJ{xB,t) 


dxf 


< 


aJ{xB,t) 


Xb 


The QFl then conforms to the inequality 


FQ{XB,t) < 


e a^J{xB,t)^ 

1 — Xd 


(Cl) 


(C2) 


The equality is obtained when the attenuation factor 
is an homogeneous function of degree a, i.e. strictly 
obeys a power law. The maximum of the QFl, re¬ 
gardless of the kind of control applied, is then obtained 

= Jo = 1 -f 

ilF(—2e“^) ~ 0.8, where W{z) is the Lambert function 
which by definition satisfies z = W {z)e^^^'’ for any com¬ 
plex number z. When the control is such that the equality 
in (C2) is satisfied, then the optimal total control-time 
at which the measurement should be done, topt, is such 
that J{xBitQpj^) — Jo- 

Under this condition, the resulting ultimate bound for 
the relative error in the estimation of Xs, which holds for 
arbitrary control on the probe, is 


when 


dJ(xB,t) 

dxB 


1 J f — e 

e{xB,t) > --- rr = — -1= 

XB \/NmFQ{xB,t) CtV Jm Jo 6 aJNm 

(C3) 

with £o = ^_w(-2e-2)(l+iW(-2e-2)) ~ ‘^■^ 8 . 


1. Precision bounds for tbe key parameters g and Tc 

(a) The probe-bath interaetion strength g 

The attenuation factor is an homogeneous function 
of degree a = 2 in the effective probe-bath interac¬ 
tion strength 5 , where G{tc,u!) = g^S{u}) with S{u}) the 
normalized spectral density of the environmental noise. 
Therefore, Eqs. (Cl) and (C2) are satisfied and the min¬ 
imal relative error in the estimation of g is obtained by 
measuring at the time topt, such that J{g,topt ) = Jo, if 
S{uj) is known. 
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(b) The correlation time Tc 


The derivative term in the QFI of Eq. (B4) depends 
on the derivative of the bath coupling-spectrum (spectral 
density) with respect to the correlation time ^, through 


dJjTct) 

Btc 


f 


dioFt{uj) 


dG{Tc,u;) 

dTc 


(C4) 


Spectral densities of the baths characterized by 
Gp(s){JctU}) = g^S{u}) [defined in the main text, Eqs. 


(7-8)] satisfy 


dS. 


/3(o) 


dTc 


< Specifically, 


dGp 


d ( Apg^Tc \ 

dTc 


dTc U + w^rfy 


(C5) 




1 - 


< 


l+TSUJ^ 

- 1 ) 

(l + rfo;^) 


1 + 


where the equality is attained when the spectral density 
function behaves as an homogeneous function of degree 
= /3 — 1, oc and 


a 


dGs 

dTc 


A 

Btc 


((s + l)rrV) 


Tc 


(C 6 ) 


since G^ is an homogeneous function of degree a = s + 1 
upon neglecting the cutoff region. 

Using Eqs. (C5-C6) to bound Eq. (C4), we find 


dJijct) 


dTc 


dujFt{uj) 


dG{Tc,ui) 


dTc 


(C7) 


< 


a 


dujFt{u})G{Tc,uj) 


aJ{Tc,t) 


Tc 


This leads to the tight bound for the Eisher informa¬ 
tion (C2) and therefore, for the minimal error in the 
estimation (C3) which holds for arbitrary control on 
the spin-probe. The ultimate bound in the precision 
(C3) is attained when the optimal total control-time at 
which the measurement should be done, topt, is such that 


Appendix D: Attainment of the nltimate bound 


1. Attainment of the ultimate precision bound 
under optimal dynamical control in the estimation 
of key parameters: 


Therefore, Eqs. (Cl) and (C2) are satisfied and the min¬ 
imal relative error in the estimation of g is attained by 
measuring at a time topt such that J ( 5 , topt) = Jn if ^(a;) 
is known. 

If S{lo) is not known, then some constraints apply for 
topt- Considering that the attenuation factor is given by 
(A 8 ), then, if the the filter function Ft{uj) is much wider 
than G{g,uj) = g‘^S{uj), it can be considered as a con¬ 
stant Ft in the integral (A 8 ). The attenuation factor fol¬ 
lowing the integration is then where we have used 

the normalization property of the spectral density S(uj). 
This limit holds when interval is shorter than the correla¬ 
tion time Tc of the environment, so that the qubit-probe 
evolves freely, yielding [15-17, 73, 74] 

FA®(a;) = ^sinc'^iujt) Ri y if t < Tc, (D1) 

and 

(D2) 

Then, the bound in the estimation error (C3) is achieved 
when 

.free .free 2 Jo _ 

Upt = provided = W — < t^. (D3) 


This condition can only be fulfilled for large enough Tc- 
To overcome this limitation, we may exploit dynamical 
control by means of frequent, stroboscopic quantum non¬ 
demolition (QND) measurements of the qubit-probe. If 
N QND unread measurements are performed during a 
total time t, the dynamics conforms to Zeno regime with 
an attenuation factor [15-17, 73, 74] 

= (D4) 

Then, the condition for attaining the bound in the esti¬ 
mation error (C3) is relaxed by 1/'/N, as 


xZeno 

^opt 



provided 


j.Zeno 

^opt 


N 



The latter condition can be attained for 


Tc- 

(D5) 


IV » 


2Jo 


<7^2- 


(D 6 ) 


In order to ensure that (A 8 ) is satisfied, suffice it that 
the product gTc be roughly estimated, by observing the 
change in the decay law from the anti-Zeno (AZE) or 
Fermi Golden Rule to the QZE regime as N increases 
[75]. 


(i) Probe-bath interaction strength g 


(a) Correlation time Tc 


The attenuation factor is an homogeneous function of 
degree a = 2 in the probe-bath interaction strength g. 


Control of the spin-probe by A^-pulse CPMG sequences 
or A^-cycles of continuous-wave (CW) driving, leads to 
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a filter function Ft{uj) that converges to a Fourier series 
described by a sum of delta functions (band narrow fil¬ 
ters) centered at the harmonics of the inverse cycle time, 
kuJa = /c G N, provided the total control time ex¬ 
ceeds the interval between the pulses t ^ Under these 
conditions, the attenuation factor (A8) becomes [23, 30] 


J{rc,t) = g^^Ft{kuJo)G{rc:kujo)- (D7) 

k^l 

In what follows we infer the conditions for attaining the 
bound from the attenuation factor for the bath spectra 
considered in the main text. 

(1) For generalized Ornstein-Uhlenbeck spectra Gp, one 
condition to achieve the bound of Eq. (C3) is that the 
filter only overlap with the power-law tail spectra. This 
is ensured when the first harmonic of the filter is already 
in this spectral region, which amounts to <C Tq- The 
attenuation factor (D7) then becomes 


Jp{Tc,t) 


^ irkNrc 

V t 


ApTc 

TrkNrc ^ ^ 


(D8) 


with Ap = ^sin{^) and cp = for CW 

(only the first harmonic k = 1 gives non-null terms) and 
cp = for CPMG (only odd k gives 


non-null terms) where C is the zeta function defined for 
Re{z) > 1 as ({z) = A. Both constants have similar 
values Cp^ ~ Cp^^^. 

Since Jp satisfies the equality in Eq. (C2), the ultimate 
bound (C3) is achieved when Jp{Tc,topt) = Jcs- This 
yields 


.opt 

^P 


= Tc 



2 2 ■ 

cpg^T^ 


(D9) 


The resulting requirement on N, in order to keep the op¬ 
timal narrowband approximation and overlap only with 
the power law tail 


^opt 




AT ■ 

ttjV tt 


(DIO) 


is then 


The ultimate bound (with a = 1) will then be achieved 
if the probe-state is measured at inde- 

P y 'c 

pendently of the control. This bound is always greater 
(worse) than the ultimate bound of Eq. (C3) for /3 > 2. 
For /3 = 2, where they are equal. 

(2) For generalized Ohmic spectra Gs, the attenuation 
factor (D7) is 


/ nkNTc 


/C=l 


Csg'^T^+^N^ 




^ irkNTc ^ * 
(D13) 


where Cg = - —(2^ ~ 1)^ ^ with kc = 
where the square brackets denote the integer part. Since 
the harmonics that contribute with non-null terms are 
those below the cutoff, then kc is defined so as to satisfy 
< ujc < and therefore -I- 1 < fee < 

t 

TtNTc ' 

The first requirement to achieve the bound (C3) is to 
avoid any overlap between the filter and the cutoff region, 
i.e. ^ Tc, with k G N. This is most likely to be sat¬ 
isfied under the narrowband approximation that allows 
to design the filter to be null at cutoff. 

The second requirement is that the filter should overlap 
with the power-law region. Therefore, one needs at least 
one harmonic below the cutoff that avoids the overlap 
with the cutoff, i.e. 


t 

■kNTc 


ttN 


> Tc- 


(D14) 


Both conditions are ensured under CW control if ^ > 
Tc, since the corresponding filter contains only one har¬ 
monic, leading to a simplified expression for the con¬ 
stant Cs = - —Then, Jg satisfies the equality in 
Eq. (C2) and the ultimate bound (C3) is achieved when 
JsiTctopt) = Jo, yielding 


J cgg^r^N^ 

V Jo 


(D15) 


To maintain this narrowband approximation, with one 
harmonic below the cutoff frequency, i.e. 


fV>maa;|- 

I cpg^T^TTP+^ j 

If the control is constrained (say, by maximum total en¬ 
ergy [26, 27, 76]), it may happen that the filter-function 
overlaps only with the Markovian region of the spectral 
density, G^{tc,uj) ~ ApTc, where it becomes a homoge¬ 
neous function of order a = I on Tc. The attenuation 
factor is then independent of the dynamical control and 
is given by the Markovian limit 

J^iTc,t)=g^Tct. (DI2) 


■j^opt FF 

— » ^ > Tc, 

TT TTiV 


(DI 6 ) 


the number of cycles N for super-Ohmic spectra (s > 1) 
should satisfy 


J'pjj.s —1 

dAguperOhm ^ 2 T superOhm ^ 1 

Csff Tc 

and for sub-Ohmic spectra (0 < s < 1) 

1 < iV™bO/^m < (D18) 


Cgg^T^ 
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which is attainable when ^ S^c- For Ohmic 

spectra (s = 1 ) 

1 < Nohm = (D19) 

implying that it must satisfy ^ gTc- 

When the last conditions cannot be satisfied, less re¬ 
strictive optimal solutions can be found by demanding 
the filter function to be null at the cutoff frequency Wc, 
and then finding the optimal control time. 


(b) The power law exponent 


To estimate the power law (PL) exponents s and /3 of 
spectral densities: oc or G(w) oc 

with a = const., we define the exponent xs = 7 that is 
common to both spectral densities: 7 = 1-I-/3 or 7 = 1 —s, 
we note that the derivative of the attenuation factor with 
respect to the exponent satisfies (cf. Eqs. (1), (4) in the 
main text) 


dJ{xB = 


97 


< 


' —00 

r*00 


dujFt{uj) 

/ CO 

dojFtioj) 

-00 


9G(7,w) 


dj 

dGP^{n,Lo) 


(D20) 


97 


97 


This expression tightly bound the QFI (B4).To attain 
the ultimate bound for the relative error of the estima¬ 
tion by the maximized QFI, one needs to apply a control 
that only probes the power-law regime of the spectral 
density. The attenuation factor under such control can 

be expressed as , where the 

dephasing time T 2 depends on the applied control. Then, 


97 




(D21) 


and the ultimate precision bound in the estimation is 
then given by Eq. (13) of the main text. The bound 
is achieved when the qubit-probe is measured at topt = 
with Ji = Ri 0.106 when the dynami¬ 

cal control can ensure this regime. For example, under 
CW or CPMG control for estimating the power law of an 
Ornstein-Uhlenbeck process, the time interval between 
the number of pulses should be smaller than the correla¬ 
tion time, <C Tc, so that N 


(c) Unattainability of the ultimate bound under free 
evolution (free-induction decay - FID) 

We here discuss the unattainability of the ultimate pre¬ 
cision bound when the spin-probe evolves freely, for the 


estimation of: 

(i) probe-bath interaction strength g (the conditions to 
attain the relevant bound are discussed in the main text 
starting with Eqs. (Dl) and ending in Eqs. (D3)), and 

(ii) the correlation time Tc, for which, as observed in 

Fig. 2 of the main text, the estimation precision be¬ 
comes worse as gxc grows. For the generalized Ornstein- 
Uhlenbeck spectra G^, as discussed in the main text, the 
free-evolution filter f/^^‘^{oj), Eq. (Dl), overlaps with 
the o; ~ 0 region and consequently the equality in Eq. 
(C5) cannot be fulfilled. For the generalized Ohmic spec¬ 
tra Gs, when t one can approximate Eq. (A 8 ) by 

its zero-order term, considering that is indepen¬ 

dent of w, since G(rc, oj) is much narrower than f/^‘^‘^{u!). 

The attenuation factor in this regime becomes indepen¬ 
dent of Tc, as in Eq. (D2), Jp'~''‘'{Tc,t) = kc 

^g^t^- In this regime there is no information concerning 
the correlation time Tc- If \g^T^ ^ Jo, then it will not be 
possible to achieve the bound under free evolution. The 
attenuation factor will depend on Tc only when t > Tc, 
and therefore \g'^J > ^ Jo, implying that the 

ultimate bound on Eq. (C3) cannot be achieved. 

The limitation of the estimation of for the spec¬ 
tra Gs using freely evolving spin-probes is that the only 
control parameter is the time t that may not simulta¬ 
neously avoid the overlap of the filter function with ujc 
and render the attenuation factor equal to JTq, i.e. usu¬ 
ally topt 4:7TnTc, n G N. Hence, the ultimate bound is 
generally not achieved under free evolution. 


(d) Pulse-error effects 


In the preceding analysis, we have considered ideal, 
i.e. perfect and stroboscopic, pulses. If non-ideal pulse- 
effects are important, one needs a model for their effect on 
the dephasing. Typically the pulses are applied with the 
same phase, as in CPMG, where flip-angle errors are com¬ 
pensated and the pure dephasing assumption is still suit¬ 
able [30, 77]. In this case, if finite-width pulse effects are 
significant, they can be considered by modifying the effec¬ 
tive filter functions [78]. Therefore, within this model for 
the error effects, our ultimate bound is valid and achiev¬ 
able in the presence of such pulse-error effects. If pulse 
imperfections generate an effective Hamiltonian that de¬ 
viates from a pure dephasing [77, 79] , then a different ap¬ 
proach to calculate the ultimate bound must be pursued. 
However, for parameter estimation, control pulses with 
good fidelity are essential. Alternatively, one may have to 
resort to more elaborate controls, nicknamed by us “selec¬ 
tive dynamical recoupling” (SDR), which we have shown 
to allow selective probing of the targeted noise source by 
factoring out uncontrollable noise sources, such as an in¬ 
trinsic T 2 decoherence due to a white-noise (Markovian) 
process or due to pulse imperfections [10, II]. 
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Appendix E: Information gain for the real-time 
adaptive estimation protocol 

The best total control and measurement time t' to be 
chosen for the next iteration of the real-time estimation 
protocol are determined by the averaged information gain 
from the currently available probability distribution of 
xb, i), assuming that experimental outcome d 

was obtained in the iteration measured at t [56]. 

The expected probability to obtain the outcome d' by 
measuring the controlled spin-probe at time t' from the 
current probability distribution oi Xb is 


Then, the best time for the next control/measurement 
tjn is defined by the value that maximizes the expected 
information 


U{tm) = maxf lj2p{d'\t\d,t) (E4) 

I d' 

/ p{xB\d', t', d, t)log{p{xB\d', t', d, t))dxB 


(E5) 


p{d'\t',d,t) = j p{d'\xB,t)p{xB\,d,t)d7 


(El) 


The total information gain of a measurement at time t' 
is given by 


d' 

where U {d\ t') is the information gain if the measurement 
at t' gives the result d'. 

The information gain of an outcome, according to in¬ 
formation theory, is measured by the entropy 

U{d\t')= / p{xB\d',t\d,t)log{p{xB\d\t\d,t))dxB■ 

m 
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